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Abstract 



The partly symmetric real Ginibre ensemble consists of matrices formed as linear combinations 
of real symmetric and real anti-symmetric Gaussian random matrices. Such matrices typically 
have both real and complex eigenvalues. For a fixed number of real eigenvalues, an earlier work 
has given the explicit form of the joint eigenvalue probability density function. We use this 
to derive a Pfaffian formula for the corresponding summed up generalized partition function. 
This Pfaffian formula allows the probability that there are exactly k eigenvalues to be written 
as a determinant with explicit entries. It can be used too to give the explicit form of the 
correlation functions, provided certain skew orthogonal polynomials are computed. This task 
is accomplished in terms of Hermite polynomials, and allows us to proceed to analyze various 
scaling limits of the correlations, including that in which the matrices are only weakly non- 
symmetric. 
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1 Introduction 



In random matrix theory the Ginibre ensembles [12] refer to Gaussian random matrices with 
either real, complex, or real quaternion entries, which are all independent. In distinction to 
ensembles of Hermitian matrices, the support of the eigenvalues in the Ginibre ensembles is a 
disk in the complex plane. Thus the eigenvalue distribution can be regarded as specifying a point 
process in a two-dimensional domain. In this paper we will study the point process associated 
with the eigenvalue distribution for random matrices which interpolate between the real Ginibre 
ensemble and the Gaussian orthogonal ensemble (GOE) of real symmetric matrices. 

The two-dimensional point process associated with the eigenvalue distribution for random 
matrices which interpolate between the complex Ginibre ensemble and the Gaussian unitary 
ensemble (GUE) of complex Hermitian matrices has been the subject of an earlier study |llj . 
Similarly, the eigenvalue distribution for random matrices interpolating between the real quater- 
nion Ginibre ensemble and the Gaussian symplectic ensemble (GSE) of real quaternion Hermi- 
tian matrices has also been analyzed as a point process [15j . The interpolating ensemble to be 
studied herein is thus the last of those naturally associated with the Ginibre ensembles to be 
considered from this viewpoint. Such studies are well motivated for their relevance to Efetov's 
theory of directional quantum chaos [6], in which a special role is played by the interpolating 
ensembles in the weak non-Hermiticity limit. 

The point processes associated with the Ginibre ensembles have the feature of being propor- 
tional to e~@ u for potentials U which are the sum of one and two body terms, the two body 
terms being logarithmic. As the logarithmic pair potential is that for two-dimensional charges, 
there is thus an analogy with the equilibrium statistical mechanics of certain two-dimensional 
one-component Coulomb systems (see e.g. [7J and references therein). This is most immediate in 
the case of the complex Ginibre ensemble, for which the eigenvalue probability density function 
(PDF) is proportional to e~^ u with (3 = 2 and 

1 N 

U = -^2\zi\ 2 - log\z k -Zj\, zj := xj +iyj. (1.1) 

1=1 l<j<k<N 

The potential is due to N unit two-dimensional charges, repelling via the logarithmic pair 
potential — log \ z — z'\, and with a smeared out disk of uniform neutralizing charge centred about 
the origin of charge density — 1/-7T, which creates the one-body harmonic potential \\z\ 2 . Note 
that the disk must have radius y/~N to neutralize the iV mobile charges. As we expect equilibrium 
Coulomb systems to be locally charge neutral (otherwise an electric field would be created, and 
the system would go out of equilibrium), the particle density should to leading order also be a 
disk centred about the origin of radius \/~N, a fact which can be checked upon exact calculation 
of the one-point correlation. We remark that recently variants on the eigenvalue problem for 
complex Ginibre matrices have been formulated [17], which have analogies with one-component 
Coulomb systems on the surface of a sphere [3] and in a hyperbolic disk p3] (in relation to the 
latter see too the work |19j on certain random complex polynomials). 

The eigenvalue PDF for matrices interpolating between the complex Ginibre ensemble and 
the Gaussian unitary ensemble (GUE) also has a Coulomb gas analogy. Such matrices can be 
written in the form H+ivA where H and A are Hermitian matrices from the Gaussian unitary en- 
semble (GUE), scaled so that the joint PDF of the elements is proportional to exp(— y^TrX 2 ), 
r = (1 — v 2 )/(l + v 2 ). The eigenvalue PDF can be computed as being proportional to [H] 

1 N 

3=1 l<j<k<N 
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Here the one-body potential can be interpreted as being due to a uniformly charged ellipse, 
charge density — r 2 ), semi-axes A and B given by [TOl 15] 

A = y/N(l + r), B = VN(1-t). (1.3) 

Again, one can anticipate that the particle density will to leading order be of this same shape, 
a fact which can be verified by exact computation of the one-point correlation |10j . 

As for the complex Ginibre ensemble, the leading eigenvalue support of the real Ginibre 
ensemble is also a disk [5]. Moreover, in the case of matrices interpolating between the real 
Ginibre ensemble and the GOE, it has been anticipated that the support will be an ellipse |18j . 
By exact calculation of the one-point function, verification of this fact is given in Section 5.1 
below. 

We begin in Section 2 by defining matrices interpolating between the real Ginibre ensemble 
and the GOE in terms of a linear combination of random real symmetric and anti-symmetric 
matrices. Knowledge of the eigenvalue PDF for the real Ginibre ensemble allows the eigenvalue 
PDF of such matrices to be computed exactly. The resulting functional form is of an identical 
structure to that for the real Ginibre ensemble, allowing in particular the corresponding gener- 
alized partition function to be written as a Pfaffian. This is considered in Section 3 and used 
to give the probability pk,N that an N x N (N even) member of the interpolating ensemble has 
exactly k real eigenvalues as a determinant of size N/2. The Pfaffian formula for the generalized 
partition function implies the fe-point correlation functions can be written as a 2k x 2k Pfaffian, 
with entries given in terms of skew orthogonal polynomials. The latter are with respect to the 
skew inner product corresponding to the entries of the Pfaffian. The main technical task is the 
computation of these skew orthogonal polynomials. In Section 4 simple expressions in terms 
of Hermite polynomials are obtained. The final section, Section 5, is concerned with various 
scaled limits of the real-real and complex-complex correlations. For fixed r the forms obtained 
are identical to the r = case (real Ginibre ensemble) except for a simple scaling of the coor- 
dinates which accounts for the change in the two-dimensional density. The case r = has been 
previously studied in [91 [22] and most comprehensively in the work of Borodin and Sinclair [2] 
(this latter work treats too the general real-complex correlations). The weakly non-symmetric 
limit of the correlations is also computed. 



2 Definition of the ensemble and the eigenvalue PDF 

Let S be an element of the Gaussian orthogonal ensemble of N x N real symmetric matri- 
ces, and thus have PDF of its independent elements proportional to e _Tr5 / 2 (equivalently, the 
diagonal elements have distribution N[0,1] while the strictly upper triangular elements have 
distribution N[0, l/y/2]). Let A be an element of the anti-symmetric Gaussian orthogonal en- 
semble of real anti-symmetric matrices which has PDF of its independent elements proportional 
to e TrA2 / 2 (each strictly upper triangular element is thus independently distributed according 
to N[0, l/y/2]). With < r < 1 and c := (1 — t)/(1 + t) define random matrices A according to 

X = ±=(S + VEA). (2.1) 

When r = and b = 1, X = S + A. In this case each element of X is independently distributed 
as a standard Gaussian N[0,1] and so A" is a member of the real Ginibre ensemble. When r = 1 
and b = 1, X = S and so A is a member of the GOE. Thus A interpolates between the real 
Ginibre ensemble and the GOE as the parameter r is varied from 1 down to 0. 
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The probability measure associated with the matrices S and A is 

{2^)- N l\- N ^ N -^l 2 e-^ s2 -^ A ^/\dA){dS). (2.2) 
From (|2,ip we compute that 

(dX) = 2 N ( N - i y 2 (^) N{N - 1)/ \Vb)- N \dS)(dA), 

and we observe too that S and A can be written in terms of X and X T . Thus we can change 
variables in (|2.2p to obtain for the PDF of the matrices X 

^, b exp(- ^_ (TtXX t -tTtX 2 )), (2.3) 

where 

Ar, b = (V~c)- N ^/ 2 (Vb) N2 (2nr N2 / 2 . (2.4) 

We seek the eigenvalue PDF corresponding to (12. 3p . 

A fundamental point is that because X is real, there is a non-zero probability that the 
eigenvalue will be real, and furthermore all complex eigenvalues must occur in complex conjugate 
pairs. Thus the eigenvalue PDF decomposes into a sum of PDFs Pk,(N-k)/2({^j} j=i,...,k', { x j ± 
iyj}j=i,...^N-k)/2'i T ib) corresponding to having k real eigenvalues {Xj} j=i,...,k and (N — k)/2 
complex conjugate pairs of eigenvalues {xj ± iyj}j=x 1 ...JN-k)/2 (for this to be non-zero k and N 
must have the same parity, a condition which will henceforth be assumed). In the case r = 
and b = 1 the probability Pk,(N-k)/2 has been computed explicitly in [13 [5] to give 

P k,(N-k)/2({^j}j=l,...,k': { x j ± Wj}j=l,...,(N-k)/2l °i 1) 

'A({Aj} i=lj ... )fc u {xj ± iyj}j=i,...,( N -k)/2) 



2 n(n+i)/4 y\? =1 r(Z/2) fc!((JV - fc)/2)! 

{N-k)/2 

xe -Ei=i^/ 2 e E J=1 fef-z?) erfc(v / 2yj) (2.5) 

where A({2j,}p=i j „. )m ) := Ilj<z( z i ~~ z j)- ^ n f ac t the eigenvalue PDF for general r and b is closely 
related to this functional form. 

In the case r = and b = 1 we read off from (|2.3p that the PDF of the elements of X is 
given by A 0i ie" TrXxT/2 . With X i-> VbX/(l - t) 1 / 2 this latter PDF becomes 

A 0il 6^/2 (1 _ r) -^/2 e -TO^/2(l-r) ; (2 6) 

while the eigenvalue PDF is obtained from (|2.5p by a simple scaling and so is equal to 



h N,2 {l _ r)- N l 2 P UN ^ k)/2 ({Vb\ 3 /{l - rs 1 2 y, , ,: 

{Vb Xj /(l - rfl 2 ± iVby 3 /(l - T) l ' 2 } ]=l _ {N _ k)/2] 0, 1). 

Now ()2.6p is a factor in (|2.3p while the remaining factor, proportional to exp ^^(fzTyTrX 2 ^ can 
immediately be written in terms of the eigenvalues of X. It follows from these considerations 
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that EES 



Pk,(N-k)/2({^j}j=l,...,k) {Xj ± iyj}j=l,...,(N-k)/2) T > b ) 

A i k {N-k)/2 

= ^(i-^- 1)/2 ex P (^(E^ + 2 £ u] .;i))/Lv ,,, 

({V&V(1 - rr % , ,: {Vbxj/il - r) 1 ^ ± ^/(l - rf 2 } j=1> ^ {N _ k)/2 ; 0, 1) 



2JV(JV+i)/4 njli T(Z/2) k\((N-k)/2) 



A({Aj}i =ll ..., fc u {xj ± iyj}j=i,...,(iv-fc)/2) 



(N-k)/2 



II erfc (Vr^ % )- (2 - 7) 



Integrating P k ,(N-k)/2 over Ay G R (j = 1, . . . , k) and (x^, yj) eKj (j = 1, . . . , (N - k)/2), 
where R 2 {( x ,v) G M 2 : y > 0} gives the probability p^jv say that a matrix of the form 
(12. ID has exactly real eigenvalues. The parameter 6 then scales out of the problem and so for 
convenience may be set equal to unity. A discussion of a systematic approach to the calculation 
of these probabilities is given in Section 3 below. The case k = N, when all eigenvalues are real, 
is special and can be considered immediately. Thus, comparing (|2.5p and (|2.7p one sees that 

/'v.ollA,}, ; v: r, 1) = (Vbf* ' l\^TT^) N{N ~ 1)/2 P N , ({^}j=l,-,N; 0, 1)- 

But we know from [5] that Pn,n\ r=o = 2~ N ^ N ~ 1 ^ 4 and so for general < r < 1 

b—l 

, 2 N-iV(AT-l)/4 
PN,N=[ T —) ■ (2.8) 



1 + T 

3 Generalized partition function and the probabilities pk,N 

The generalized partition function associated with the PDF (I2.7P is defined by 



k 

/oo /*oo ™ /* ^ 

dAi--- / dA fc TTu(A/) / dxidyx--- / dxr N _ k y 2 dy(N-k)/2 
-oo J-oo ;=1 JR^ ,/ R 2 

(N-k)/2 

X ]J U ( X ^y) P fc,A r -fc({Aj}j=l,...,fc;{^ ± %-}j=l,..,(AT-fe)/2;T,^)- (3- 1 ) 

In view of the functional form (|2,7p this is structurally identical to the case r = and 6 = 1, 
when a Pfaffian formula is known [21J. Consequently ()3.ip too has a Pfaffian evaluation. 

Proposition 1. Zei {pz-i(x)}; = i i ... i tv 6e a se£ of monic polynomials of the indexed degree, and 
let 



a jtk [u}= dxu(x) dyu(y)e b{x +y )/2 p j - 1 (x)p k _ 1 {y)sgn(y - x) (3.2) 

J —oo J —oo 

(3 j>k [v\ = 2i dxdyv(x,y)e b ^- x ^evic^ 



26 

-y 



x (pj-i(x + iy)p k -i{x - iy) - Pk-i(x + iy)pj-i(x - iy)) . (3.3) 
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For k, N even we have 

z HN - k)/2 [uM = ^^^^ [C fc/2 ]Pf ItmM + Pj,Mj,i=K~,N, (3.4) 

where [C P ]/(C) denotes the coefficient of £ p in /(C)- 

We remark that the Pfaffian operation applies to even-dimensional anti-symmetric matrices; 
the formula f|3.4[) therefore requires modification for N odd. Such a modification is known |21j . 
but to avoid having to consider separately the cases ./V even and N odd, only the case N even 
will be considered hereforth (it is planned to address the case N odd in a separate publication). 

From the definitions, ^jt (jv-fc)/2[lj 1] = Pk,N and so we have 

< / I -l T )N(N-l)/2 

PkN= T, v [C fc/2 ]Pf [C«i I [1] + A J [1]] i 1=1 N (3.5) 

2 N ( N+l )/ 4 Yl^_ 1 T(l/2) '"'' 6=i 

(it is valid to set 6 = 1 since as noted below (|2.7|) Pfc,7V is independent of h). Suppose for 
definiteness that in (|3.2[) we choose Pj(x) = x J = 0, . . . ,N — 1). Changing variables x \— * —x, 
y ^ —y shows 

«2j,2fe[l] = a2j-i,2fc-i[l] = 0, (3.6) 
while we see by introducing polar coordinates and changing variables 6 ir — that furthermore 

&i,2fc[l] =/Vi,2*-i[l] = 0. (3-7) 

The equations (|3.6p and ()3.7[) give that the matrix in (|3.5p has a checkerboard pattern of zeros. 
Recalling the general formula (Pi A) 2 = det^4, rearranging rows and columns in A so that all 
non-zero entries are in the top right and bottom right N x N blocks, and using the fact that the 
entries are anti-symmetric shows the Pfaffian can be written as a determinant of half its size. 
Thus 

( / l + T ) N (N-l)/2 

Pk ' N = 2*Wi)/4 Ui i W2) ^ det[C " 2j - 1 ' 2fc[1] + ^i,2fc[%fc=i,...,*/2 (3-8) 

where the quantities in the determinant have 6 = 1. It remains to evaluate these quantities. 
In relation to OL2j—l,2k W\ > integration by parts shows 

a 2 i-i,2fc[l] = 2{k - l)a 2i _i )2fe _ 2 [l] + 2F(j + k- 3/2) 

and thus we obtain the explicit formula 

^•-i, 2fc [l] = 2*(fc - 1)! X: r ^^~ • M 

p=i ^ 

For /32f_i2fc[l] we see from the definition that 

/%-i,2fe[l] = " 4Im ^ + dxdyey 2 - x \rfc(^^ y yx + iy) 2 ^\x - iy) 2k 

2j-22fc-l 9 -_ 9 9j l_-i 

= " 4 E E (^) (^T 2 ) ("W + * - 1 - (' + P)/2K«* (3-10) 



2fc-l 



/=0 p=0 
l+p odd 
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where 



ij := J y leiic {\[^~~) ev2 d y (i ° dd )- 

Integration by parts shows 



J 2 3 J^M 1-tVI + t/ 2 u/ 1 2 J ' 



This recurrence has solution 



which makes /?2j-i,2fc[l] explicit. In the case r = (|3.1ip reduces to a result of Edelman [5]. 

In the case r = explicit formulas for pfc jv were first calculated in [5] , by direct integration 
of (|2.5p • A formula closely related to (|3.5p was given in [16] , while the formula (|3.5p itself in the 
case r = was given in [9]. A determinant formula equivalent to (|3.8p in the case r = 0, but 
derived using a different strategy, is given in [1]. 



4 Correlations and skew orthogonal polynomials 

As realizations of X are not conditioned on the number of real eigenvalues k, the generalized 
partition function Zn[u, v] appropriate for calculation of the correlation functions is obtained 
by summing (|3.ip over all allowed k, 

N 

Z N [u,v]= Z kj ( N _ k y 2 [u,v\. 



k=0 
k even 



It follows from (BHD that 



L ' J 2 Af ( iV + 1 )/ 4 r(//2) 

Correlation functions can be calculated from this by functional differentiation. For example, the 
n-point correlation function between real eigenvalues at x%, . . . ,x n is given by 

1 5 n 

p\ n) {xi,...,x n ) = — — -Z N [u,l] (4.2) 

Z N [1,1\ du(xi) ■ ■ ■ 5u(x n ) -=1 



u=l 



Moreover, all ra-point correlation functions can be expressed in terms of a 2k x 2k Pfaffian 
[2] (see also [9] in the case of all real eigenvalues in the correlation, or all complex eigenvalues). 
The structure of these formulas is the same for all allowed r and b, which is a consequence of 
the structure of the entries of (|4.ip being the same for all allowed r and b. In particular, in the 
case (IP) 

-I r (xj,Xk) S r (xj,Xf : ) 
-S T (x k ,Xj) D T (xj,x k ) 



p T {n) (xi,...,x n ) = Pf 

where with 



(4.3) 

**(*):=/ sgn(x -y) Pk (y)e-y 2 / 2 ^dy (4.4) 
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one has 



U k 



y 



k=0 

d ~ r 1 rv 

D ( x >y) = ~5Z S { x iV)i I T (x,y) = -sgn(y - x) - S T (x,z)dz. 

J X 



dx~ ' y ~ ia ' 2 

We note too the explicit Pfaffian formula for the correlation between complex eigenvalues 



(4.5) 
(4.6) 



p\ n) {{x uyx \ (x n ,y n )) = ]J (2.e^ 2 -^)/( 1+ ^erfc( J r 

l=i V 

S°(zj,z k ) S^(zj,z k ) 
S^(zj,z k ) S c T {zj,z k ) 



rVi 



xPf 



where Zj := Xj + iyj and with q 2 j_ 2 (z) := -p 2 j-i(z), q 2 j-i{z) :=p 2 j- 2 (z), 



N 



S c T (w,z)=J2 



Pj-i( w )lj-ii. z ) 
u [(i-i)/2] 



(4.7) 



(4.8) 



The dependence on r comes in through the requirements of the polynomials {Pj(z)}j=o 
In addition to being monic of the appropriate degree as in Proposition [IJ they must be skew 
orthogonal with respect to the skew inner product associated with the matrix in (|4.ip for u = 
v = l. 

Explicitly, this inner product reads 



with 



(/,9)r := 



dx 



dye 



if, 9) ■= (f,9)r+(f,9)c 

-0(z 2 +y 2 )/2f( x )g(y) sg n(y - x) 



(4.9) 



(/, g)c ■= %i J^ 2 dxdy e b{y2 x ^ eTic^^^-y) (j{x + iy)g{x - iy) - g(x + iy)f(x - iyj) . 

The set {Pj(sc)}j=o,l,... is sa id to be skew orthogonal if 

(P2j,P2k) = (P2j+i,P2k+i) = (j, k = 0, 1, . . . ) (p 2j ,p 2k+1 ) = (j,k = 0,1, ... j ^ k), 

while (p 2 j,p 2 j+i) = Uj 0. The main technical task then is to compute these polynomials. 
Note that the parameter b acts as a scale of the coordinates, and so there is no loss of generality 
in setting b to a specific value. It turns out that a convenient choice is b = 1/(1 + r). Making 
this choice, the skew inner product of interest reads 



(f,g) ■= (f,9)r + (f,g)< 



(4.10) 



with 



(f,9)r 

(f,g) c ■-- 



dx 



2i 



dx 



o 



dye f(x)g(y)sga(y - x), 

y 2 -x 2 ( I 2 \ 

dye i+T erfc I J - _ y 1 [f(x + iy)g{x - iy) - gix + iy)f[x - iy)\ 
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Theorem 1. Introduce the scaled monic Hermite polynomials 



Cn{z) = {^) n ' 2 H n {^=). (4.11) 



~2J "Ky/fr. 
The family of monic polynomials {Rj(z)}j = o i ... with 

R2n+i(z) = C 2n+l {z) - 2nC 2n -i(z), R 2n (z) = C 2n (z) (4.12) 

are skew orthogonal with respect to the skew inner product 10[ ). Furthermore, the normaliza- 
tion r n is given by 

r n := (R 2n , R 2n+l ) = (2n)!2^(l + r). (4.13) 

Proof. Note that R 2n+ \(z) is an odd polynomial while R 2n (z) is even. These properties are 
sufficient for the derivation of (13.61) and (13.71) so it is immediate that 



(R2j,R2k) = (R 2j+ i,R2k+i) = 0. (4.14) 

It remains to verify that 

(R 2j+1 ,R 2k ) =0 (4.15) 

for k 7^ j, and that for k = j the normalization (|4.13l) results. This will be done by computing 
the explicit form of the skew inner product between C 2 j+i and C 2k , 

(C w ,C^ ) = \-^ k+m ^{ k + ^^ (4.10) 
,0, j < k. 

Assuming (|4.16p . the explicit formulas (|4.12p show that (|4.15p is valid and furthermore give the 
normalization (|4.13p . The task is thus reduced to proving (|4.16p . For this, repeated use will be 
made of the properties of the Hermite polynomials (|4.1ip 

^C n (z)=nC n - 1 (z), (4.17) 

zC n (z) = C n+1 (z) + n<_i( Z ). (4.18) 

Let us first define 

roo [-co x 2 +y 2 

Ij,k= dx dye ^)C 2j+1 (x)C 2k (y)sgn{y-x). (4.19) 

J — oo J— oo 
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Then a partial integration over x gives 



/•OO y 2 

/ dye 2 ^C 2k (y) 

J GO 

j /"CO /" 

+ —rL dX L 



2j + 2 B vy ; 



dye -^T ^±gM c 2fc ( y ) S gn( y - x) 



+ f° dx /°° dye-^%±^C 2fc (y)25(y - x) 

f°° x 2 +y 2 

dx / dye 2 ( 1+T ) C 2j+3 (x)C 2k (y)sgn(y - x 



1 1 

" 1 • t 2. / ■ 2 , 

+ —— / dx / dye 2(i+r )C ' 2i+1 ( :r )c 2fc ( 2/ ) S g n (y_ x ) 

^ ~r ^~ ./ — oo J — oo 

2 r 00 z 2 

+ 2jT2 /_ ^ 1+rC Wz)C 2fc (x) 
~ l + r2j + 2 J ' +1,fc + 1 + r j ' k + 2j + 2 1 + r' 

where 

/OO ^2 
dxe" ~ C 2i+2 (a;)C 2fc (x) . 
-oo 

Similarly a partial integration over y gives 



dxe ^i+^C^+^x) 



C 2(l + r) 



2A; + 1 



sgn(y - x) 



y=co 



y=-oo 



+ - 



1 



1 + T 

coo 



<■!•'■/ ^.li-rj ^-'^ sgni// x) 

1 C , 2fc+i(y) 



dye 2 ( 1 +'")C2j+i(a 



2fc + 1 



2<% - x) 



/OO PC 
dx 
-oo J — c 
1 1 T 00 /*°° * 2 +y 2 

l + T2k + 1 J dx J d V e 2(1+ ^C 2j+1 (x)C 2fe+2 (y)sgn(y 

r f 00 //°° _£W 

/ dx / dye 2 (i+-) C 2i+ i(x)C 2fc (y)sgn(y - x) 

J — oo J — c 



1 + r 



OO </ — OO 

2 ~ 2 



2k + 1 

1 I 



dxe i+-C 2j+ i(x)C 2fc+ i(x) 



l + r2fc + l /j ' fc+1+ l + r /j ' 



k — 



where 



Vj,k 



2k + ll + T 

dxe i+-C 2j+ i(x)C 2fe+ i(x). 

-OO 



Thus we obtain recursion relations 



Ij+i,k = (2j + 2)Ij,k ~ 2^,fc, 
Ij,k+i = {Ik + l)I j:k + 2rj jik . 
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Let us next derive the recursion relations for 

/•oo rcc y 2_ x 2 

J jjk = 2i dx dye i+- erfc(7y) 

J-oo JO 

x [C 2 j+i(x + iy)C 2k (x - iy) - C 2 k(x + iy)C 2 j+i(x - iy)} (4.25) 

with 



7=^3—^. (4-26) 
For that purpose, we consider an integral 

/oo rco 2 2 

/ y — x 

dx / dye X + T erfc(7y) 
oo Jo 

x [C 2 j+2(x + iy)C 2 fe-i(x - iy) - C 2k -i(x + iy)C 2j+2 (x - iy)] . (4.27) 
For k > 1, a partial integration over x gives 

/•OO ^2 

= 2i dye 1 + T erfc(7y) 
■/ o 



/ , ■ ^C2fc(x - iy) 
e 1 +TC 2 j+2{x + iy)- 



2k 

C 2k (x - iy) 



1 f°° f 00 y 2 - X 2 

+4i / dx dye 1 + T evic( , yy)xC 2 j+ 2 (x + iy) 

1 + T J-oo Jo 

-2i j°° dx J™ dye^eMiy)(2j + 2)C 2j+1 (x + iy) ^ 
+c.c. 

2i 1 Z" 00 Z" 00 a 2 -* 2 
= — — — / dx dye erfc(-yy)xC 2j+2 (x + iy)C 2k (x -iy) 

K l + T J_ OQ J 

j + 1 f°° f°° y 2 -* 2 

-2i— — / dx dye i+- eric(-fy)C 2j+1 (x + iy)C 2k (x - iy) 

& J-oo Jo 

+c.c. (4.28) 
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On the other hand, a partial integration over y gives 
K jtk 

= 2i I dxe - " 



OO 



ei+^erfc(7y)C 2j+2 (x + iy) — 



1 r r ...... 4^....r,.,. ,.. , ,.„,< g»(x-iy) 

2k 

C2k(x - iy) 



1 r 00 r 00 y 2 -^ 2 
4i— — / dx / dye ^ erfc^y^C^+^x + iy)- 
l + r j.qo Jo 



-oo 

(■OO 



2k 



+2i / dx dye i-h- erfc(jy)(2j + 2)C 2 j+i(x + iy) 

J -oo JO 

+2^ fdx dy^^C 2]+2 {x + iy) ^F + c.c. 

V 71 " ./-OO JO 2fe 

+- / dxe C 2 j+2{x)C 2 k(x) 
ft J -co 

2i 1 f°° f°° y 2 -* 2 
— — / dx / dye i+r eiic(~{y)iyC 2 j+ 2 (x + iy)C 2k (x - iy) 
K l + T J_ 00 J 

j + 1 /■=» /-OO ^2 

+2i— ^— / dx / dye erfc(7y)C 2j+ i(x + iy)C 2k (x - iy) 

J — oo J 

2 / 2 f°° f°° x 2 v 2 

~fcV 7r(l - r 2 ) ./ dx y dye 1+T % ~ T C2 3+^( x + ^J^kix - iy) 

+c.c. (4.29) 



Comparing (I4.28P and f|4.29|> . we obtain 
2t 1 



/■OO r-oo y 2_ x 2 

/ dx / dye eifc(jy)(x + iy)C 2j+2 {x + iy)C 2k {x - iy) 
J-oo Jo 

j _|_ l r°° f°° y 2 -x 2 

-4i— 7— / dx / dye~^erfc(7y)C2 j+ i(x + iy)C2fe(x-iy) 2fc 

+ C.C. 

2 f 00 * 2 
= - / dxe ^ C 2j+2 (x)C 2k (x) 

ft J-oo 



fc 1 + T 



dx / dye X + T !- t 

-oo JO 



feV tt(1-t 2 ) 

x [C 2 j+2(x + iy)C 2k {x - iy) + C 2i+2 (x - iy)C 2k (x + iy)] . (4.30) 

Therefore, noting (|4.18p and the orthogonality relation (see e.g. [10J) 

f°° f°° y J_ . 

/ dx / dye i+^e i-t C m (x + iy)C n (x - iy) = nml V 1 - r 2 5 mi „, (4.31) 

J —oo J —oo 

we can derive 

J j+ l,k = (2j + 2) Jj- k + 2& * - 2^(2j + 2)! (1 + r) S j+1 , k . (4.32) 
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In order to derive another recursion relation, we similarly employ partial integrations to find 

roc 

= 2i / dye 1 + T erfc(7y) 
Jo 



and 



e 1 + T ^ L —. : C 2k+1 (x-iy) 



2j + l 

+4i— — / dx / dye^+r evfc^y) 2j+1 . xC 2k +i(x - iy) 

l + T J-oo JO 2 J + 1 



-2i 



dx j dye 2 fe ! erfc( 7 y)^±i^±^(2fc + l)C 2k (x - iy) 
lo 47 + 1 



+c.c. 

2i + 

-2i / dx / dye ^ erfc(7y)C 2 j+i(x + iy)C 2k (x - iy) + c.c. 

47 + 1 J-oo Jo 



1 f°° f°° y 2 -x 2 

-— — / dx dye i+- exic(^y)C 2j+ i{x + iy)xC 2k+1 {x - iy) 

1 1 + T J-oo JO 



(4.33) 



j-i,fe+i 

/•°° re 2 

2i / dxe !+ T 

J — oo 

1 



■ jil f , ,C 2j+ i(a; + iy)^ , . ,1 y= °° 
ei+^erfc(7j/j — ,. n , . — C 2fc+ i(x - zy) 



z(2j + 1) 



j/=0 



-4t 



oo poo 

t dx 

1 + 7" J-oo JO 

oo 



x 2 Co ■ _i_ i fx ~\~ iy) 

dye i+T erfc(7y) ^ - ^ yC 2fc+ i(x - iy) 



+2i 



/°° dx /°° dye4+f ! erfc( 7 y) ° 2j+ ^ X ± ^ ( 2 fc + l)C 2fe (x - /y) 

J-oo JO 



2j + l 



27 



_ (7 j /)2 C 2j+ i(x + iy) 



+2i—= I dx dye X + T e 

vW-oo 7o *(2j + l) 

2 Z" 00 z a 
- 9 . - , / dxe i+-C 2j+ i(x)C 2 fe + i(x) 

47 + 1 J-oo 

4i 



C 2fc+ i(x - iy) + c.c. 



1 f 00 f 00 ^-q; 2 

-— — / dx / dye i+- erfc(jy)C 2j+ i(x + iy)iyC 2k+ i(x - iy) 

1 1 + T J-oo JO 
2k + 1 f°° f°° y 2 -* 2 

+2i——— dx dye evfc(jy)C 2j+l {x + iy)C 2k (x - iy) 

47 + 1 J-oo JO 

4 / 2 /Z 00 /Z 00 * 2 y 2 

+ 2j + 1 V tt(1-t 2 ) J dx Jo dyG 1+T 1 ~ tC ^ +1 ^ x + i v) c 2k+i{x - iy) 



+C.C. 



(4.34) 
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A comparison of (|4,33p and (|4.34j) yields 
Ai 1 



/oo roo 2 2 

I y — x 

dx / dye X + T erfc(7y) 
-oo Jo 



2 j + 1 1 + T 

xC 2j+1 (x + iy)(x - iy)C 2k+1 (x - iy) 

2k + 1 f°° f°° y 2 -* 2 
-4i—— dx dye i+- erfc(7y)C 2j+ i(x + iy)C 2k (x - iy) + c.c. 
Z J + 1 J-oo Jo 

4 f°° x 2 
= - „. -. / dxe ^C 2j+1 (x)C 2k+1 (x) 

4 / 2 r°° , r 00 , 

H W — ; or / dx / dye 1 + T 1 ~ T 

2 j + iy vr (i- r 2)y_ oo y 

x [C 2j+ i(x + iy)C 2k+ i(x - iy) + C 2 j+i(x - zy)C 2fc+ i(x + iy)] . (4.35) 
As before it follows that 

Jj,k+i = (2fc + 1) J ijfc - 2T3P,-, fc + 2>/M2j + 1)! (1 + r) ^ fc . (4.36) 
Let us employ the notation 

= + Jj,k- (4.37) 
Then, from (|4.24j) . (j4.32[) and (14.36p . we obtain the recursion relations 



L j+ i, k = (2j + 2)L jtk - 2V2vr(2j + 2)! (1 + r) 5 j+1>k , j > 0, fc > 1, 

L jtk+1 = {2k + + 2^(2j + 1)! (1 + r) S jtk , j>0, k > 0. (4.38) 

Using these recursion relations and noting 

L , = -2^(1 + r), (4.39) 

we can readily find that (|4.16p holds. □ 
We remark that the two key Hermite polynomial properties f)4. 1T|) and (|4. 18j) allow us to 
verify that the first of the two formulas in (|4. 12[) can be rewritten to read 

i? 2n+1 ^) = -(l + r)e z2 / 2 (^)^(e- 2 / 2 ( 1 ^)C7 2 „(z)). (4.40) 

Use will be made of this form in the analysis of the correlations (|4.3p . 

Taking the limit r — * in Theorem [1] gives the family of skew orthogonal relevant to the real 
Ginibre ensemble, a result which was announced and made use of in [9]. 

Corollary 2. The family of monic polynomials {pj(z)}j=o,i,... specified by 

P2n+i(z) = z 2n+1 - 2nz 2n -\ p 2n (z) = z 2n (4.41) 

are skew orthogonal with respect to the skew inner product ^4-10 ) in the case r = 0. The 
corresponding normalization is given by 

u n :=(p 2n ,p 2n+1 ) = (2n)\2V2^. (4.42) 
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5 Asymptotic properties of the correlations 



5.1 Eigenvalue support 

The eigenvalue support for the real Ginibre ensemble (r = case) is to leading order a circle of 
radius y/N ' . To gain some insight into its expected form for < r < 1, consider the portion of 
(|2.7|) which is dependent on {zj := Xj+iyj}j =1 ^ N _ k y 2 - For yj large this portion is proportional 
to (jl.2p with N i— > (N — k)/2. As remarked below the latter equation, previous analysis of the 
one-point correlation has revealed that for the PDF (|1,2|) the density is supported on an ellipse 
with semi-axes A and B given by (jl.3p . The exact results obtained above can be combined with 
the analysis of [TO] to verify that this result persists in the present setting. 

In [10] the boundary of the support is characterized by the values of (x, y) which maximize 
the difference 

P(l){(x,y))\N^N+l - p c {1] ({x,y)) (5.1) 
for large N. We know from (f4"7?]) and P~8j) that 



p c {1) ((x,y)) = 2ie^ 2 - * 2 )/(i+-) erfc ( yj -^y) ^ (z, . 



(5.2) 



To compute the asymptotic form of (|5,ip . following the beginnings of a strategy used to analyze 
the two-point correlation for (jl.2p in [8], use will be made of an integral form of (|5.2p . In this 
regard, from a standard integral representation of the Hermite polynomials we have 

i poo 

C n (z) = — e- t2 (z + V2^it) n dt. (5.3) 

V n J-oo 

It follows from this and Corollary [5] that 

2 f'OO ^ yOO ^ 

S c T {w,z) = — -/ dhe-*! dt 2 e- t *£%(w + \fcit u z+ y/5?it 2 ). (5.4) 

7T(I + T)j-oo J-oo 

On the other hand, substituting (|4.41|) in (|4.8|) gives for Sq the simple expression 

Substituting (|5.5p in (|5.4p and making further use of (|5.3p it follows that for large N, x, y 

pli)((x,y))\N^N+i - pl 1} ({x,y)) 
V2 1 



cxp 



vr(l+r) (JV-2)! 

This same function of x, y and results from studying the difference (|5.ip in the case of the 



PDF (jl.2p . As remarked above, working in [10J deduces from this that the boundary of the 
support is given by an ellipse with semi-axes specified by (jl.3p . 

5.2 Density of real eigenvalues 

Next asymptotic properties of the density of real eigenvalues will be considered. According to 
d and (S31) 



-x 2 /2(1+t) N / 2 ~ l 1 

P\i)(x) = 2^(1 + T ) J^(®2k{x)R 2k+1 (x) - <5>2k+i(x)R2 k (x)y (5.6) 
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The mean number of real eigenvalues is obtained by integrating p T ,^ (x) over the real line. Making 
use of (|5.6p . (|4.4p (with p k replaced by R k ), (|4.40|) and Theorem [1] shows 

Further, use of a tabulated integral |X3|. §7.373], and a Kummer transformation for 2-P1, gives 

I" e- 2 ™ 2 /( 1+ ^(F 2fc (x)) 2 dx = 2 2fc - 1 / 2 (i±^) 1/ V 2fe - 1 /2r(2A ; + 1/2) 

x 2 Fi(l/2, 1/2; -2k + 1/2; -r/(l - r)). (5.8) 

For large the %Fi function is to leading order equal to unity. Hence the leading order behaviour 
of a general term in (|5.7p is 

/l + r\i/2 1 



V^rd -r) 



(2fc)V2 



and so for large N 



To analyze the density itself for large N, we again make use of (|4.4|) (with replaced by 
as well as the identity 

R2k + 2(x) - (2k + l)R 2k (x) = -(1 +^) el2/2(1+T) ^( e " l2/2(1+T) ^ + i(s 
(cf. (|4.4U|) ; this can be verified using (|4.17|) . (|4.18|) ) to rewrite (|5.6|) in the simplified form 

= ~W g A?! ^ + 2^(1 + r) (iV-2)! ' (5 - 10) 
Use can now be made of the classical summation formula 

tkH k{x)H k {y) i 2j-l/2 e -*a(x»-Hf , )/(l-t a ) e 2xwt/(l-ta) j | t | < j (5 n ) 

fc=0 

to conclude 

pffi ( x ) ■= lim ^ n (x) = - 1 (5.12) 



With the leading order support of the real eigenvalues the interval [— v N(l + r), v JV(1 + r)], 
to leading order the mean number of eigenvalues must be equal to 2y/N(l + r) (the length of 
this interval) times the density (|5,12p . This reclaims (|5,9p . 

We turn our attention now to the neighbourhood of the spectrum edge. In the case of the 
real Ginibre ensemble (r = 0) the explicit form of the density profile about the spectrum edge 
at x = y/~N was exhibited as [9] 

Jta/fafc^ + JOL = 4-(i(l - «f V/2X) + ^(1 + erf X)). (5.13) 
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For general < r < 1 the density at the spectrum edge is analyzed by setting x = (l+r)y/N +X 
in (|5,10p then taking the limit N —* oo. As is distinct from the bulk scaling, the summation and 
the term distinct from the summation both give 0(1) contributions. Consider first the latter. 

To calculate the explicit form of the contributions, our main tool is the Plancherel-Rotach 
asymptotic formula [2UJ, [3] 

rr r \ m \n ( %2 - aV^ 2 - 2n - n r^. — \ 

H n [x) ~ (2n) exp y nlog(x — y ar — 2n)J 



2 \ \/x 2 — 2n 



valid for n large and x > \/2n. Recalling ()4.1ip we require this formula with x = yN/2(y r r + 
1/y/r) + X/\^2t, n = N — k. A straightforward but tedious calculation gives that for k fixed 

H N - k (x) ~ (2N) N - k e- k (l-rr 1/2 

/ _/Y / \ X 2 k k \ 

x exp [t - log 2N - log rj + VNX - — _ - + k + - log 2N + - log rj .(5.15) 

Consider now the final term in (|5. 10j) . The asymptotic form of C^-i(x) can be read off (|5. 14j) 
by setting k = 1. To use it to deduce the asymptotic form of &n-2(x) we require the integral 
evaluation [13] 

e- x2 H 2m (xy) dx = 0F^(y 2 - 1) 



e ri2m\xy) ux = v vi — yy - i v " ! 

-OO 

This formula allows us write 



N-2 



/ t nAT/2-1/ (AT -2)! , v/ , 



e-'"/2P+T) ffiv _ 2 ( » j rft j. (5. 16 ) 



It is in this form that we substitute (|5.15p with k = 2. Combining results and making use too 
of Stirling's formula allows us to compute the sought limiting form, 



lim 



JV-oo V2v^F(1 + t) {N-2)\ 

To analyze the sum in (|5.10p the asymptotic expansion (|5.15p must be extended to include 
terms 0(k/y r N). One finds these terms to be the multiplicative factor 

k 2 T kX 

exp 



2JV(1-t) VN{1-t)S 
Noting too from Stirling's formula that 

' (2vrAf) 1 / 2 exp (A^logA^ - N- klogN + — V 



(N — k)\ y J "V & " " & ' 2NJ- 
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after rearranging the order of summation so that k *— > N — k, {k = 2, . . . , N) and recognizing 
that a Riemann sum approximation to a definite integral results we find 



-*V(1+t) ^ 2 i 
lim (C k (x)y 



x=(l+r)y/N+X 



Now adding together (|5.17p and (|5.18p gives for the edge density 



== i/faCP + rK» + x) - V2j(1 _ t2) 

x(i(l - erf(V2X/(l - T 2 ) 1/2 )) + e " A 2 ^" r ' > (l + erf(X/(l - r 2 )" 2 ))) .(5.19) 

Note that this agrees with (|5.13|) in the case r = 0. 

We observe from (j5.19j) that p^ e (X)dX for general < r < 1 is obtained from the case 

t = by the simple scaling X Xj\J\ — r 2 . Note that this same rule is valid for the bulk 
density of real eigenvalues (j5. 12|) . Indeed it is reasonable to expect that all local correlations are 
only altered by this change of scale, as the point process for general < t < 1 is locally identical 
to that for the point process in the case r = 0, except that the two-dimensional bulk density is 
scaled by a factor of 1/(1 — r 2 ). We will now proceed to exhibit this fact for the general real-real 
and complex-complex correlations in the bulk. 

5.3 /c-point correlations in the bulk 

Consider first the complex-complex case. Setting 

S c T (w,z) := e~^ +w2 ^ 2( - 1+ ^S^(w,z), 
we see that (|4.7p can be rewritten 



I 2 

n(^erfc(J T — ^^Pf 
i=i v 



S°(zj,z k ) S^(zj,z k ) 



(5.20) 



Now, it is immediate from (|5.5p that 



lim Sq(w,z) 



N-+oc ' 2^/2^ 
Substituting this in f)5.3j) and computing the resulting Gaussian integrals gives 



lim S c T (w,z)= }lL ^ expf- . T 9N {z 2 + w 2 ) + —L-^zw). (5.21) 

N^oo tV ' 2^(1 -r 2 ) V 2(l-r 2 ) V ' I - t 2 ) V ' 

Consequently 

^ *<»• 2) - i^rS) exp ( - wS) ■ (5 - 22) 
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Substituting in (|5,20p gives the bulk complex-complex correlations. The feature that the bulk 
limiting value of 

n 

P(n)(( x i>yi)'- • • d x n,y n )) Y[dxidyi (5.23) 

1=1 

for general < r < 1 is gotten from the r = case by the replacements 



(xi,yi) i ^ (xi/ \J\-t 2 , yi/\/ 1 - r 2 ) 

is evident. 

It remains to consider the real-real case. Proceeding as in the derivation of (|5,10p shows that 
(14.51) can be rewritten 

e -(xW)/2(l+r) y 2 /2(1+r) 9 f x ) 

s '<- »> ^ £ ^<^-« + 2 vs (1 + T) (5 ' 24) 

As for the derivation of (|5.12p . we use (|5.1ip to both deduce that the final term vanishes as 
N —* oo (a consequence of the convergence of the sum) and to give a closed form evaluation of 
the summation. It follows that 

lim S r (x,y) = -= 1 e -(x-y) 2 /2(i-r=') (5 _ 25) 

In view of the formulas (|4.6|) for the remaining quantities in (|4.3|) , as for (|5.23p we have that the 
limiting bulk value of 

p T ( n) (xi, . . . , x n ) dxi ■ ■ ■ dx n 
for general < r < 1 is gotten from the r = case by the replacements Xj i— ► Xj/\Jl — r 2 . 

5.4 The weakly non-symmetric limit 

In relation to the ensemble interpolating between the complex Ginibre ensemble and the GUE, 
it was exhibited in [11] that well defined correlations result by setting r = 1 — a 2 /N, then taking 
N — > oo. Similarly, scaled correlations were computed in this limit for the ensemble interpolating 
between the real quaternion Ginibre ensemble and the GSE. Here the scaled correlations of the 
real Ginibre/ GOE interpolating ensemble will be calculated. Note that with this scaling (II. 3p 
gives that the eigenvalue support collapses onto the interval [—2y/N, 2y/N] of the real axis. The 
mean spacing between eigenvalues is then 0(1/ \/~N), suggesting that we should also multiply 
coordinates by tt/^/N (the proportionality ir is chosen for convenience; then a unit real density 
results) before taking N — > oo. 

Now, using the asymptotic expansion 

^fi±Ae~ x2 H n (x) = cos ( v / 2nTTx - mr/2) + O^ 1 ' 2 ), 
1 (n + 1) V / 

we deduce from (|5.24p that 

7T / ttx iry \ 1 [F^e-^l:! / IT, \ 
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Making use of Stirling's formula, a Riemann sum approximation to a definite integral is obtained, 
and we compute 



lim * S*'™ ^ 



N VN 



T =l-a 2 /N 



cosiru(x — y) du. 



(5.26) 



Note that in contrast to the correlations implied by (|5.25|) . we see from (|5.26p that the correla- 
tions in the present setting of the weakly non-symmetric limit exhibit an algebraic decay. For 
the limiting form of the complex-complex correlations, we first note from (|5.3p - ([5.5p that 



S C T (W,Z) 



2(1 + r)V2vr 



N-2 

£ 

3=0 



C j+1 (w)Cj(z) - Cj(w)C j+1 (z) 

r(i + 1) 



Proceeding now as for the working which lead to (|5.26p shows 



lim , 



TT 



si 



TTW TTZ 



N VN 



T =l-a 2 /N 



TT 

— I ue 
2 



siinru(w — z) du, 



(5.27) 



where on the LHS we have used the fact that the complex-complex correlations must be scaled 
by (tt/^/N) 2 for each independent two-dimensional coordinate (x,y) to account for the measure 
in (|5.23p . Substituting this in (|4.7p and noting too that 



,(«)/(l+T) erfc 



1-T 2 



III 



<rl ' C( v) 



gives the explicit weakly non-symmetric limiting form of p c , v 

[n) 
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